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The spin and charge excitations of the t-J ladder are studied by exact diagonalization techniques for 
several electron densities. The various modes are classified according to their spin (singlet or triplet 
excitations) and their parity under a reflection with respect to the symmetry axis along the chains 
and a finite size scaling of the related gaps is performed. At low doping formation of hole pairs 
leads to a spin gap for all I/t ratios. This phase is characterized by (at least) one vanishing energy 
mode only in the charge bonding channel when q x — > consistent with dominant superconducting 
correlations. At larger doping the spin gap disappears. Although the anti-bonding modes remain 
gapped, low energy q x ~ and q x ~ 2kF spin and charge bonding modes are found consistent with 
a single band Luttinger scenario. At sufficient low electron density and above a critical value of J/t 
we also expect another phase of electron pairs with gapped spin excitations. 
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The physical properties of t-J ladders are of inter- 
est both because of their possible realization in various 
materials such as (VO)2P2C>7 1 and SrCu203 2 and be- 
cause they offer new insights into the nature of magnetic, 
charge density, and pairing correlations in strongly cou- 
pled electron systems. Here we report results for the 
charge and spin excitation spectrum of a doped isotropic 
t-J ladder obtained from exact diagonalization calcula- 
tions on 2 x 10 and 2x5 clusters. Based on these results, 
we discuss a possible J/t-^n^) ground state phase dia- 
gram. 

The t-J Hamiltonian on the 2 x L ladder is defined as, 

H = J' ^JS . • S i>2 - \n jtl n jt2 ) 
i 

P,j 

-tJ2 P G(c},i3, s Cj+i,P,s + H.c.) J P G 

-t'^^(c] )1;s c J)2;s + H.c.)P G , (1) 

j> 

where most notations are standard. [3 (=1,2) labels the 
two legs of the ladder (oriented along the «-axis) while 
j is a rung index (j=l,...,L). In the following we shall 
consider the isotropic case where the intra-ladder (along 
x) couplings J and t are equal to the the inter-ladder 
(along y) couplings J' and t' respectively. However, to 
clarify the discussion, we have explicitely separated the 
intra- and inter-ladder parts in the Hamiltonian. 

The half-filled case (Heisenberg ladder) is well 
understood. 3-6 When J=0, % reduces to isolated rungs 
and the spin excitations can be constructed by exciting 
singlet rungs to triplets each costing an excitation en- 
ergy, J'. Once the intra-ladder coupling J on the rungs 
is turned on, the triplets can propagate, form a coherent 



band (with a minimum energy at momentum K x =ir) 5 , 
and the spin gap is reduced. 

Upon lightly doping the t-J ladder the formation of 
hole pairs with persistence of the spin gap has been estab- 
lished by various numerical 7,8 or analytical techniques 9 . 
The nature of pairing is particularly evident when J' is 
the largest energy scale since, in this case, it is energeti- 
cally favorable for each hole pair to lie on a single rung 
to minimize the energy cost due to the destruction of 
singlet bonds. Interestingly enough pairing continues to 
exist in the isotropic case J = J', even at small ratio 
J/t or in the regime of intermediate couplings described 
by the Hubbard model 10 . These bound pairs were also 
found to be quite robust under increasing the number 
of coupled chains. 11 However, the evolution of the spin 
gap under further doping is still unclear. At low electron 
density a phase of electron pairs exists above a critical 
ratio J/t where the two-electron bound state becomes 
stable. Therefore, under hole doping, a transition from 
a hole pair phase to an electron pair phase must occur. 
Since at small electron or hole concentrations the spin 
gaps of the paired phases are continuously suppressed by 
electron or hole doping respectively there is a possibility 
that an intermediate phase with gapless spin excitations 
may exist. 

The goal of this paper is then to identify the various 
possible phases by determining the exact nature of their 
low energy excitations. A semi-quantitative derivation 
of the phase diagram can then be attempted. The low 
energy modes are first characterized by their spin. Sin- 
glet and triplet excitations correspond to charge and spin 
modes respectively. It is also very usefull to consider the 
parity of the excited states under a reflection with re- 
spect to the symmetry axis of the ladder along the di- 
rection of the chains. Even (R x = 1) or odd (R x = — 1) 
excitations correspond to bonding or anti-bonding modes 
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respectively. 12 Lastly, the total momentum K x along the 
ladder is considered in order to determine the form of the 
dispersion relation of each mode. 

Before going further it is intructive to consider the ex- 
citation spectra for free fermions on a ladder. 13 At low 
enough doping the fermi level intersects the anti-bonding 
and bonding branches (split by an energy 2t' = 2t) at 
±^1 and ikp t 2 respectively. Excitations of bonding 
(anti-bonding) nature correspond to intra-band (inter- 
band) transitions. The bonding excitation spectrum con- 
sists then of four gapless modes (both in the spin and 
charge sectors) at momenta q x ~ (2 modes), 2kp t \ 
and 2hp t 2- In the anti-bonding sector vanishing energy 
modes appear only at finite momenta ±(^i ± kp^)- It 
is already clear that this simple picture breaks down in 
the t-J ladder in the region where the spin excitations 
are gapped. When the hole pair concentration is very 
small one can clearly identify two different kinds of spin 
excitations; 8 (i) a triplet excitation of a singlet rung away 
from the hole pairs (typical of the half-filled case) with 
R x = — 1 and (ii) an excitation of a hole pair into 2 prop- 
agating quasi-particles with R x = 1. Crudely, these two 
processes involve a finite excitation energy proportional 
to J'. Note that, with increasing doping, these bonding 
and anti-bonding spin modes loose their simple physical 
interpretation. At intermediate electron densities when 
the Fermi level only crosses the anti-bonding branch (for 
(n) < 0.5) we shall see that the weakly interacting pic- 
ture acquires nevertheless some relevance. In that case 
the anti-bonding (interband) modes are all gapped and, 
as far as the low energy modes are concerned, one only 
deals with a single band. 

As seen above, numerical techniques have already 
proven to be very successful in studying coupled chains. 
Exact diagonalizations (ED) methods are particularly 
well adapted to the investigation of the low energy modes 
since implementation of quantum numbers and symme- 
tries is rather straightforward. The following study is 
based on exact diagonalisations of 2 x 5 and 2x10 dou- 
ble rings at intermediate electron densities (^n^ = 0.4 and 
(^n^ = 0.8 corresponding to the less understood region of 
the phase diagram. Note that the electron number is al- 
ways N e = Ap to guaranty that the antiferromagnetic cor- 
relations are not frustrated when one goes around each 
chain. Special care is also needed when choosing the 
boundary conditions (BC). Finite size corrections are re- 
duced by using BC such that a non-interacting Fermi sea 
(obtained by turning off the interaction) would form a 
closed shell. This is achieved with anti-periodic (peri- 
odic) BC for (n) = 0.4 ((«) = 0.8). The various ex- 
citation modes are then straightforwardly obtained by 
calculating the GS energy in each symmetry sector with 
the Lanczos algorithm. 

Charge and spin modes are shown in Figs, l(a-d) for 
(ri*l = 0.8 as a function of J/t. Data for all possible 
momenta are displayed but the precise behavior of the 
excitation energies vs K x will be discussed later on. Let 



us rather concentrate first on the energy gaps defined 
as the smallest excitation energies for each set of quan- 
tum numbers (ie the 4 possible bonding and anti-bonding 
charge and spin excitations) irrespective of momentum. 
In Figs. 1 the lowest excitation energy (gap) of the 2 x 10 
cluster at each J/t value is indicated by a full square. 
The long-dashed line passes through the equivalent gap 
for the 2x5 cluster at the same electron density. Under 
doping the system becomes metallic and general hydro- 
dynamic arguments predict long wavelength zero energy 
charge fluctuations. We then expect the finite size gap 
associated with the bonding charge mode to scale like 
1/L (critical behavior) 14 . The short-dashed line in Fig. 
1 corresponds to a simple a + bj L extrapolation of the 
2x5 and 2 x 10 data (where it is not shown the extrap- 
olation gave a small negative number). As seen in Fig. 
1(a) this procedure was found to give a very small (or 
even negative) gap in agreement with the existence of a 
zero energy mode in the R x = 1 (bonding) charge chan- 
nel. When we examine the K x momentum dependence 
we will see that this corresponds to the long wavelength 
K x —> mode. We expect that this type of extrapo- 
lation will in general give a lower bound for the gap 14 
since, in the true asymptotic regime (ie when the system 
size exceeds the correlation length), the finite size correc- 
tions change to an exponantial behavior with 1/L. Figs. 
1(b) and (d) strongly suggests the existence of a gap in 
the anti-bonding spin and charge modes. Fig. 1(c) shows 
that a gap also exists above a small critical value (our 
data give J/t\ c ~ 0.4 but J/t\ c could be much smaller). 
We also find that all the gaps scale almost linearly with 
J/t. Above J/t\ c , the relevant physical picture is then 
consistent with a gas of hole pairs propagating along the 
chain direction. Below J/t\ c it is likely that the hole pairs 
dissociate into two quasi-particles (QP) of opposite spins. 
In this phase a zero energy bonding spin excitation can be 
simply obtained by flipping the spin of one of the quasi- 
particles. However, an anti-bonding excitation requires 
the transfer of an anti-bonding QP to a bonding QP or 
the excitation of a singlet rung into a triplet. Both of 
these processes require a finite energy cost in agreement 
with Figs. 1(b) and (d). 

Hole pairs are also expected to become unstable with 
increasing doping. Similar data as Fig. 1 for the charge 
and the spin excitations are shown in Figs. 2(a-d) for a 
smaller electron density (n) = 0.4. Again Fig. 2(a) is 
consistent with the existence of zero energy (long wave- 
length) charge fluctuations characteristic of a metallic 
behavior. Figs. 2(b) and (d) show the existence of a siz- 
able gap in the spectra of the anti-bonding modes. This 
is physically clear in the picture of the non- (or weakly) 
interacting fermi gas where the Fermi level only intersects 
the lower anti-bonding band (for (n) < 0.5). Electron- 
hole transitions from the lower band to the upper band 
are then gapped and we expect the gaps to depend weakly 
on J in agreement with Figs. 2(b) and (d). The interpre- 
tation of Fig. 2(c) is more subjective. The extrapolated 
gap remains small and, contrary to Fig. 1(c), is weakly J 
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dependent. This suggests that the spin gap vanishes in 
this region of the phase diagram. 

More insights on the nature of the two spin gap and 
gapless phases found previously can be obtained from 
the dispersions of the four modes (ie the excitation en- 
ergies vs K x ). The latter are shown in Fig. 3(a) and (b) 
for (n) = 0.8, J/t = 2 and (n) = 0.4, J/t = 0.5 re- 
spectively. In Fig. 3(a) the low frequency charge fluctu- 
ations are clearly separated from the three other modes 
at higher energy. Note also that there is no indication 
for any finite momentum zero mode contrary to the non- 
interacting case. Nevertheless, an investigation of the dy- 
namical charge correlation function would clearly be use- 
ful to study the number of R x = 1, K x —> charge modes 
(other zero bonding modes with larger charge velocities 
could also exists). Fig. 3(b) is more intricate. The bond- 
ing charge and spin excitation branches are very similar 
and seem to exhibit q x ~ and q x ~ 2k p low energy 
modes while both anti-bonding modes are gapped. This 
is reminiscent of a single band Luttinger Liquid 15 . 

We conclude this paper with a brief discussion of a pos- 
sible phase diagram which is sketched in Fig. 4. At large 
J/t the system separates into a hole-free and a hole-rich 
phase. The phase separation line at large J/t has been 
computed by Tsunetsugu et al. 8 from the compressibil- 
ity. At smaller J/t ratios, when the electron density is 
reduced our data suggest a transition from a spin gap 
phase to a Luttinger-like phase 15 with gapless spin and 
charge excitations. At small electron density another 
paired phase must also appear due to the existence of 
a two-particle boundstate above J/t\ c = 2. Other 2p- 
particle (p > 1) boundstates could also become stable at 
larger J/t in this region. 

In the spin gap phase, we find evidence that there is 
only one gapless charge density mode which occurs for 
K x —> in the R x = 1 (bonding) channel. This is con- 
sitent with the power law behavior of the singlet rung- 
rung pairing correlations found in recent numerical renor- 
malization group calculations (NRG). 10 Furthermore, the 
gap in the charge excitation spectrum at 2k p is also con- 
sistent with the NRG results which fail to find the conju- 
gate, Luther-Emery, 16 power law decay of the 2k p charge 
density correlations. It would appear that these conju- 
gate correlations must developp in another channel be- 
sides the particle-hole charge density excitations studied 
here. 
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FIG. 1. Charge ((a) and (b)) and spin ((c) and (d)) ex- 
citations of a 2 x 10 ladder at electron filling (n^ = 0.8 
as a function of the ratio J/t. The bonding (anti-bonding) 
modes correspond to the figures (a) and (c) ((b) and (d)). 
The energy reference is given by the GS energy and periodic 
boundary conditions in the chain direction are used. Data 
for all possible momenta (K x = 2irn/L, n integer) are shown 
(open and full squares) and the various data sets correspond- 
ing to the same momentum are connected by a continuous 
thin line. The lowest excitation energy (gap) is indicated 
by full squares for the 2 x 10 cluster and is also shown for 
a smaller 2x5 cluster (at the same electron density) by a 
long-dashed line. The short-dashed line corresponds to a ten- 
tative 1/i-extrapolation of this finite size gap. 
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FIG. 2. Same as Fig. 1 for an electron density of = 0.4 
and with anti-periodic boundary conditions in the direction 
of the chains. 

FIG. 3. Spin and charge excitation modes of a 2 x 10 ladder 
versus momentum q x (in unit of 7r). The quantum numbers 
associated to the various symbols are shown on the plots, (a) 
and (b) correspond to (n^ = 0.8 and J/t = 2 and to (n^ = 0.4 
and J/t = 0.5 respectively. 

FIG. 4. Tentative schematic picture of the J/t - (n^ phase 
diagram. The boundary with the phase separated region is 
taken from Ref. 8 



